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Abstract. Identifying dark matter and characterizing its distribution in the inner region of 
halos embedding galaxies are inter-related problems of broad importance. We devise a new 
procedure of determining dark matter distribution in halos. We first make a self-consistent 
bivariate statistical match of stellar mass and velocity dispersion with halo mass as demon- 
strated here for the first time. Then, selecting early-type galaxy-halo systems we perform 
Jeans dynamical modeling with the aid of observed statistical properties of stellar mass pro- 
files and velocity dispersion profiles. Dark matter density profiles derived specifically using 
Sloan Digital Sky Survey galaxies and halos from up-to-date cosmological dissipationless 
simulations deviate significantly from the dissipationless profie of Navarro-Frenk- White or 
Einasto in terms of inner density slope and/or concentration. From these dark matter pro- 
files we find that dark matter density is enhanced in the inner region of most early-type 
galactic halos providing an independent dynamical evidence for halo contraction. The main 
characteristics of halo contraction are: (1) the mean dark matter density within the effective 
radius has increased by a factor from ~ 1 for clusters with M^ir > lO^^M© to 4 - 5 for 
galaxies with Mvir < IO^^Mq where Mvir is the halo virial mass, and (2) the enhancement is 
more frequently realized by steepened density slope than increased concentration compared 
with the fiducial NEW profile. Based on our results we predict that halos of nearby ellip- 
tical and lenticular galaxies can be promising targets for 7-ray emission from dark matter 
annihilation. 
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1 Introduction 

The currently standard cold dark matter (to be referred to as ACDM in conjunction with 
Einstein's cosmological constant A) model of structure formation requires that galaxies form 
at the centers of dark matter halos [1, 2]. Dark matter halos embedding galaxies are indicated 
by rotation curves of spiral galaxies [3-5], velocity dispersions of early-type (i.e. elliptical 
and lenticular) galaxies [6-8], and gravitational lensing effects induced by galaxies [9-11] 
and galaxy clusters [12, 13]. In the ACDM model initial small density fluctuations grow by 
gravity and overdense regions collapse to form halos dominated by CDM and halos merge 
subsequently to form larger and larger halos in a hierarchical fashion. These halos provide 
gravitational potential wells for baryonic matter to sink dissipationally to form visible galaxies 
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in their centers. Current ACDM-based galaxy formation models reproduce well observed 
clustering properties of galaxy distributions and are ever improving towards the goal of 
explaining all observed statistical and intrinsic properties of galaxies [14-16] . 

Despite the fact that dark matter halos are pivotal for the standard (ACDM) model of 
structure formation to explain observed galaxies, our current knowledge of them is limited in 
two main aspects: dark matter particles have never been identified conclusively [17-19] and 
statistical properties of dark matter distribution within halos have not been characterized 
fully and robustly through observational studies [20-23], hydrodynamic simulations [24-28] 
or empirical modeling of galaxy formation [29-31]. In particular, detailed dark matter den- 
sity profiles including both inner density slope and concentration have been obtained only 
for a limited number of individual systems [32-37]. Precise knowledge of the dark mat- 
ter distribution within halos is crucial for the search of dark matter particles because both 
direct and indirect detections of dark matter particles depend on the dark matter density 
within halos [19, 38]. The dark matter distribution within halos is expected to be modified 
by dissipational galaxy formation in their centers from the initial distribution predicted by 
dissipationless N-body simulations [24-27, 39, 40]. Hence, a precise statistical characteri- 
zation of dark matter distribution is crucial for galaxy formation physics as well [24, 25]. 
Galaxy formation modeling requires properly taking into account halo contraction effects. 
This means that observationally determined contracted halo mass profiles provide useful 
constraints on galaxy formation physics including dissipational gas cooling/accretion, star 
formation/feedback, supernovae and AGN. 

Here we present a new procedure of determining the distribution of dark matter within 
the halos embedding early-type spheroidal galaxies. We intend to produce a rigorous and 
general statistical characterization of halo mass profiles that is allowed by currently best 
statistical knowledge of galaxies and halos. Our approach is to combine a statistically repre- 
sentative sample of observed galaxies with halos from a cosmological dissipationless N-body 
simulation. Specifically, we have two complete and separate sets of data for galaxies and 
halos: photometric and spectroscopic data for galaxies [41, 42] from the Sloan Digital Sky 
Survey (SDSS) [43] and halos from the Bolshoi simulation [44] (section 2). Based on the data 
sets we follow two steps to infer the dark matter distribution in halos. The first step is to 
make a statistically consistent one-to-one match between galaxies and halos for early-type 
systems by simultaneously assigning stellar masses (M^) and velocity dispersions (a) to halo 
virial masses (Mvir) so that observed statistical functions (of and a) and correlations 
(between and a and between and Mvir) are reproduced (section 3). Each galaxy with 
and a assigned is further assigned the effective radius (Re) and the Sersic index (n) of 
the light distribution using their observed correlations with and a (appendix A). The 
next step is to carry out Jeans dynamical modeling of each system with the assigned galaxy 
parameters and the unknown dark matter distribution (section 4) . Based on the stellar mass 
profile and velocity dispersion we constrain only the allowed range of dark matter density 
profiles for each system (section 4.3.1). Out of the allowed range we determine the most likely 
dark matter density profile for each halo using the statistical property of currently available 
velocity dispersion profiles within (section 4.3.2), which is then tested against physically 
and observationally well-motivated dynamical constraints (appendix C). We discuss the im- 
plications of our results for halo contraction (section 5) and dark matter annihilation strength 
(section 6). We conclude in section 7. 
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2 Data: galaxies and dark halos 



According to the ACDM paradigm a galaxy is embedded near the center of its host halo. 
The halo itself may be embedded in a larger halo or may embed smaller halos referred to as 
subhalos. The larger host halo and the subhalos also embed galaxies at their centers. In this 
way one can identify a one-to-one match between galaxies (including satellites) and halos 
(including subhalos) [15, 45-50]. A reliable match requires unbiased samples of galaxies and 
halos. 

2.1 Observed galaxies from the SDSS 

A large and unbiased galaxy sample is provided by the completed SDSS. This survey covers a 
large area of sky providing a low redshift galaxy sample unaffected by large scale structures of 
the Universe [43]. Galaxies are observed not only photometrically but also spectroscopically 
providing structural and kinematical parameters such as the total stellar mass M^, the stellar 
velocity dispersion u, the effective radius i?e, [41] and the Sersic index n [42]. The veloc- 
ity dispersion a refers to the luminosity-weighted line-of-sight velocity dispersion (LOSVD) 
within Rc/8 as described in [41]. Velocity dispersions are actually measured within a fixed 
aperture of radius 1.5 arcsec. The fixed aperture corresponds to different physical scales that 
are smaller than for most galaxies. Measured velocity dispersions are then corrected to a 
fixed physical radius of i?c/8 using a typically observed radial profile of luminosity- weighted 
LOSVDs within i?e- This corrected velocity dispersion thus corresponds to the luminosity- 
weighted LOSVD within Re/8. 

Galaxies can be classified by eye inspection, color, spectral features, or luminosity profile. 
None of these methods are perfect in separating galaxies by morphology. With the goal of 
dynamical modeling we select a galaxy population that consist mostly of early-type spheroidal 
systems which are predominantly velocity dispersion supported systems. Our selection is 
based on a luminosity profile concentration index Cr, which is the ratio of the radial scale 
which contains 90 percent of the Petrosian luminosity in the r-band to that which contains 
50 percent. Spheroidal (early- type) galaxies are known to have more concentrated profiles 
compared with disk (late-type) galaxies. We use the criterion of Cr > 2.86 to select spheroidal 
galaxies [41]. Our working assumption is that for our selected galaxies any rotating disks, if 
present, can be ignored for the Jeans dynamical modeling. 

2.2 Theoretical halos from the Bolshoi dark matter only simulation 

Cosmological N-body simulations can be used to produce samples of dark matter halos. The 
predicted statistics of these halos without galaxies are robust classical results of cosmological 
physics. The predicted statistical properties of halos, however, depend on the adopted values 
of cosmological parameters. We use the Bolshoi simulation [44] that is based on the following 
up-to-date cosmological parameters of the flat ACDM cosmology: h = 0.7, Om = 0.27, 
Hs = 0.95 and as = 0.82 consistent with Wilkinson Microwave Anisotropy Probe seven-year 
results [51]. The Bolshoi simulation differs from the Millennium simulation [14] particularly 
in that the latter adopts a higher erg = 0.9. The Bolshoi simulation uses 2048'^ particles in a 
cosmological box of size 250 Mpc allowing a statistically representative complete catalog 
of halos with virial mass (defined below) 1.5 x lO^°/i"^M0 < M^ir < (1 - 2) x lO^^/i'^M©. 

The Bolshoi simulation identifies both isolated halos (distinct halos) and satellite halos 
(subhalos) embedded in larger halos and characterizes their properties separately. For any 
mass range the halo population is dominated by distinct halos (from more than 70 percent 
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to 100 percent) but we include subhalos for accuracy. The Bolshoi simulation provides the 
distribution of halos in maximum circular velocity. The maximum circular velocity is found 
to be tightly correlated with the virial mass Myir that is defined to be the mass bounded by 
the virial radius within which the mean density is equal to the virial overdensity Avir times 
the mean cosmic density = ilm/Ocrit (/Ocrit being the critical density of the Universe). For 
the above flat ACDM cosmological model Avir ~ 360 and the virial radius rvir is related to 
Mvir by 



at redshift z = Q [52] . The halo mass function then follows from the combination of distinct 
halos and subhalos at fixed Mvir- The Bolshoi simulation does not sample halos beyond 
Mvir ~ (1 — 2) X 10^^/i~^Mq. However, an analytic extrapolation of the halo mass function 
beyond the mass limit agrees with a larger volume simulation ( 'MultiDark') results [53] 
based on the same cosmological parameters. 

3 Statistically matching galaxies with halos: a bivariate log-normal distri- 
bution of stellar mass and velocity dispersion as a function of halo mass 
for early-type systems 

Abundance matching of a galaxy parameter (usually luminosity or stellar mass) with halo 
mass has been widely used recently [15, 45-50]. For the early-type galaxy population we carry 
out a rigorous abundance matching of two parameters (stellar mass and velocity dispersion) 
of a galaxy with its host halo mass. There will result a bivariate distribution of stellar mass 
and velocity dispersion as a function of halo virial mass. The halo virial mass includes all 
mass within the halo, i.e. Mvir = + where is the dark matter mass and is 
the stellar mass of the galaxy.^ The radial dark matter density profile within the halo with 
M(jm is the unknown that we want to solve for. The galaxy is characterized primarily by its 
stellar mass and velocity dispersion a. The radial stellar mass density profile within the 
galaxy with is modeled by a deprojected form of the Sersic mass profile [54]. The effective 
radius Rc and the Sersic index n can be assigned according to the observed correlations with 
and a (appendix A). Thus, the primary task is to assign simultaneously and a to 
Mvir for all halos in a statistically consistent manner. 

The successful match between galaxies and halos (i.e. the assignment of M* and a to 
Mvir) requires that the results satisfy all known statistical functions and correlations. These 
include the halo mass function, the galaxy stellar mass function, the galaxy velocity dispersion 
function and the observed correlations between Mvir and M^ and between M^, and a? The 
simultaneous assignment of and a to Mvir can be done through a bivariate probability 
distribution of M^ and a as a function of Mvir. The observed distribution of cr at fixed (or 
vice versa) is well described by a log-normal distribution for the early-type galaxy population 
[41, 55]. The observed distribution of M^ at fixed Mvir can also be described by a log- normal 
distribution [56]. These observations justify a choice of a bivariate normal distribution for 
the logarithmic values of M^ and a at fixed Mvir- Note that a bivariate normal distribution 
is not an accurate model for the total galaxy population, which includes both spheroidal 

^ We ignore the mass of interstellar gas that is non-negligible for spiral systems because our work is concerned 
only with spheroidal systems. Black holes harboring at the galactic centers are also negligible for our analysis 
although we include them for completeness. 

^Note that no direct observational constraint is available for the correlation between Afvir and a. 
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and disk galaxies, because the observed distribution of a at fixed for all galaxies is not 
log- normal but asymmetric. 

We derive the bivariate normal distribution for the logarithmic values of and o" as a 
function of Mvir for the early-type galaxy population using an iterative procedure described 
in detail below. The parameters of the bivariate normal distribution are the mean Myir-M^ 
relation, the mean Mvir-o" relation, the standard deviations of log^gC-^*) ^-^d log;^o('^)) ^^'^ 
the correlation coefficient between logj^gl-^*) log^Q((T) as functions of Mvir. The mean 
relations are based on an abundance-matching relation between Mvir and M^ and the SDSS 
observed relation between M^ and a. Starting with initial guesses of the standard deviations 
and the correlation coefficient we iterate keeping the mean relations fixed until the bivariate 
distribution reproduces the SDSS galaxy stellar mass and velocity dispersion functions and 
the M^-o" relation up to the observational uncertainties. 

3.1 The detailed procedure 

We describe in detail the procedure^ of generating a bivariate log-normal distribution of stellar 
mass (M*) and stellar velocity dispersion (a) as a function of halo virial mass (Mvir), or a 
bivariate normal distribution of y = logiQ{Mi, / Mq) and Z = logiQ(c7/km s~^) as a function 
of X = logio(^vir/MQ). We use a sufficiently large comoving volume of 4 x 10^ Mpc^ to keep 
high precision in statistical match toward the large mass limit. Our results are then limited 
only by the accuracies of the input statistical quantities. Throughout a system refers to the 
combination of a halo (whether it is isolated /distinct or a subhalo embedded in a larger halo) 
and the central galaxy embedded at the center of the halo. 

The statistical quantities of our use for galaxies are (1) the stellar mass function (SMF) 
for all-type galaxies, (2) the SMF and the velocity dispersion function (VDF) for early-type 
galaxies, (3) the number fraction of early- type galaxies Je as a function of M^, and (4) the 
distribution of early-type galaxies in the stellar mass- velocity dispersion plane, to be referred 
to as the M^-cr relation. The M^-a relation for early-type galaxies with Cr > 2.86 is derived 
and presented here while all the other have been published [41]. Stellar masses are based on 
the Chabrier initial mass function throughout. The statistical quantity of our use for halos is 
the halo mass function (HMF) that includes both distinct halos and subhalos. Throughout 
a statistical function (i.e. SMF, VDF, or HMF) is denoted by (p and expressed as comoving 
number density per unit logarithmic interval of the variable under consideration. 

1. The first step is to find a mean relation between and Mvir for all- type galaxies 
from the SMF and the HMF using the abundance matching method [15, 45-50]. This method 
uses a prior knowledge (or assumption) of the intrinsic scatter of one variable at the other. 
Observational studies find sy = 0.16-0.17 for the intrinsic scatter of Y at fixed X insensitive 
to the value of X [56-58]. We adopt sy = 0.17 [56]. We follow the procedure demonstrated in 
[50] . Initially we obtain an approximate abundance matching relation assuming zero intrinsic 
scatter. This initial relation is biased at large X owing to the ignored scatter. The difference 
between the initial relation and the unknown true relation is to be referred to as bias. The 
remaining task is to estimate the bias and thus the true mean relation as well. We do 
this iteratively using a Monte-Carlo method in the following way. Let us first consider the 
initial mean relation as a surrogate of the true mean relation. We generate halos from the 
HMF and then assign Y to X through the surrogate mean X-Y relation and the intrinsic 
scatter sy = 0.17 assuming the Gaussian model. From these mock systems we derive a mock 

^The reader who is only interested in the result may skip this lengthy subsection. 
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Figure 1. The first panel sliows the halo mass function from the Bolshoi N-body simulation. The 
solid curve includes all halos (i.e. both distinct and subhalos) while the dashed curve includes only 
distinct halos. Data points show mock halos within a comoving volume of 4 x 10^ Mpc'^. The second 
panel shows the SDSS galaxy stellar mass function for all-type galaxies. The curve is the SDSS 
measured function while data points are from galaxies assigned to the halos according to the relation 
shown in the last panel. The last panel shows the final abundance matching M^i-^-Mi, relation with a 
Gaussian intrinsic (constant) scatter of 0.17 for log(A/*) at fixed Mvir- 



SMF which is of course quite different from the observed SMF because of the biased input 
surrogate relation. Now obtain a mock abundance matching relation between this mock SMF 
and the HMF ignoring the intrinsic scatter. Then the difference between this mock relation 
and the input surrogate relation is our first estimate of the bias. Using this bias we correct 
the initial abundance matching relation between the observe SMF and the HMF. With the 
corrected relation we do the simulation all over again. There results a better corrected 
relation. We iterate this simulation until the mock SMF matches the observed SMF. Usually 
a few iterations suffice. Fig. 1 shows the reproduced HMF and SMF and the final X-Y plane. 
The total number of the generated systems with Mvir ^ lO"^"'^ Mq is ~ 1.866 x 10^. 

2. We use the above X-Y plane to convert the observed early-type fraction of systems as 
a function of Y (stellar mass) to that at fixed X (halo mass). Let /e(-'^, Y) be the early- type 
fraction as a function of X and Y. If /e(X, y) were known, /e(X), the fraction at X would 
be simply / f^{X,Y)P{Y\X)dY where P{Y\X) is the probability distribution of F at X 
that comes from the X-Y relation. In our case P{Y\X) is the Gaussian distribution with a 
standard deviation of sy = 0.17. Without fE{X,Y) we use the following approximation 
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Figure 2. (Top) A schematic view of converting the early- type fraction as a function of Y to that 
at fixed X. We use /e(X) = E^[MX,Yi)w, + /e(X, « EJ/e(^.) + /E(r-z)]w;, where u., 

and w^i are statistical weights satisfying + W-i) = 1- The last approximation follows from 

u)_i = and the assumption that {[/e(X, - fE{Yi)] + [/e(X, - /E(F_i)]}w^i ~ 0> which is 
plausible (see the texts). (Bottom) The left panel shows the early- type number fraction (black curve) 
of SDSS galaxies as a function of stellar mass. We use the observed fraction of galaxies with Cr > 2.86 
for lO^'^M© < M^, < since outside this range the measured fraction is less reliable [41]. The 

red galaxy number fraction [41] is also shown for comparison. The right panel shows the early-type 
number fraction (black solid curve) as a function of halo mass (Mvir) converted from the observed 
fraction at Mi, using the Mvir-M^ plane shown above. The dotted curve is a linear extrapolation as 
a function of log]^o(Mvir) for the high-mass part (Mvir > lO^^'^M©). The gray curve is the number 
fraction of visually selected early- type galaxies based on halo occupation statistics results [42] . The 
red curve shows the converted red number fraction for our red galaxies. The hatched regions show 
observational results (95% confidence regions) on the red fraction based on SDSS data and satellite 
kinematics (green/blue region based on the observed fraction as a function of luminosity/stellar mass) 
[56]. There are consistencies for both early-type fractions and red fractions. 

where /e(^) is the observed fraction at Y. This approximation is obtained as follows: 
/e(X) = I fE{X,Y)P{Y\X)dY 

^ Y,[fEiX,Yi)P{Yi\X) + MX,Y^i)P{Y^i\X)]AY 

i 

= Y,[h{X,Yi)wi+_f^iX,Y_i)w.ii (3.2) 




loglo(Mvir/Mo) 

Figure 3. Halo mass function for those halos embedding early- type galaxies at their centers (red 
curve) derived by multiplying the total halo mass function (black curve) by the early-type fraction 
(shown in Fig. 2) as a function of Mvir- 



where Yi [i = 1,2,3, •••) are greater than the mean (y), y_j = 2(Y) — Yi (see the top 
panel of Fig. 2), and statistical weights Wi[= P{Yi\X)AY] and W-i[= P{Y^i\X)AY] satisfy 
^i{wi + W-i) = 1 and W-i = Wi as well for the symmetric Gaussian distribution. It then 
follows 

h{X) = Y,[h{X,Yi) + h{X,Y.i)]wi 

i 

= Y.^H^^) + h{y-^)]wi + Y}.^h{X,Yi) + 5h{X,Y_i)]wi, (3.3) 

i i 

where fE{Y) = f f{X,Y)P{X\Y)dX [P{X\Y) being the probability distribution of X at 
fixed Y] and 5hiX,Y) = h{X,Y) - fE{Y). But, we expect 5h{X,Yi) « -5h{X,Y.i) 
for Yi close to iy) (i.e. small values of i) assuming that f-E{X,Y) is a smooth distribution 
and the shape of P{X\Y) is slowly varying with Y . For large values of i the approximation 
5f-E[X, li) ~ — (5/e(^, Y-i) may not be good enough but w for the Gaussian distribution 
of Y . Hence, we have X^J'5/e(-'^; ^i) + 5f'E{Xy-i)\wi and equation (3.3) is reduced to 
equation (3.1). 

Fig. 2 shows the functional behaviors of observed and converted fractions. The early- 
type [Cr > 2.86) fraction is compared with the red fraction. At a given stellar mass the 
red fraction is significantly higher than the early-type fraction. This implies that significant 
fractions of red galaxies do not possess concentrated light distributions. The converted 
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Figure 4. Same as Fig. 1 but for early-type galaxy-halo systems. 



fractions as functions of Mvir are compared with recent independent results. Notice that 
our early-type fraction based on luminosity concentration index is in good agreement with 
a fraction based on visual inspection of galaxy morphology [42]. Our red fraction is also 
consistent with independent results based on color selection [56]. We obtain the HMF for 
early- type systems by multiplying the number density for all- type halos by /e(X). Fig. 3 
shows the HMF for early- type systems. 

3. Assign y to X for the early-type population using the same procedure of step 1 
based on the early-type HMF from step 2 and the observed early-type SMF. Fig. 4 shows 
the resulting HMF and SMF and the X-Y plane for the early-type population. We obtain 
~ 6.9 X 10^ early-type systems with X and Y assigned. 

4. For the early-type systems from step 3, we assign Z to Y using only the observed 
Y-Z relation (Fig. 5) ignoring the correlation between X and Z. In this way we obtain an 
approximate set {X,Y, Z}. This set is approximate in the sense that parameter correlations 
and scatters have not yet been determined. Nevertheless the mean relations among three 
parameters should be accurate as long as the distributions in the two-parameter planes are 
assumed symmetric [observations support the symmetric distributions of Y (at fixed X) and 
Z (at fixed Y) for the early-type population]. This set allows us to obtain a crude estimate 
of sz (standard deviation of Z) at fixed X. 

5. Let us now assume that Y and Z follow a bivariate normal distribution at fixed X. 
The probability density function is given by 



P{Y, Z\X) = ^ exp 

2-KSYSzJl - p\z 



2(1 - P\z) 



(3.4) 



with 

[y-PY? , [Z-^iz? 2pYz{Y - py){Z - pz) 
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Figure 5. (Left) The black curves are the SDSS steUar mass {M^,) and stellar velocity dispersion (a) 
functions for early-type galaxies with Cr > 2.86 [41]. Red points are derived from the data points 
shown in the right panel. (Right) The black full and dashed curves represent our measured median 
values and standard deviations for SDSS early-type galaxies based on the data analyzed in [41]. Red 
data points show a sample of mock galaxies produced from our final bivariate distribution of and 
(7 as a function of Mvir to be shown in section 3.2. The purple filled circles and upright/upside-down 
triangles are respectively the median values and standard deviations of the data points. 



where /iy and fiz are the mean values and pYZ is the correlation coefficient at fixed X. 
The values of fiy and fiz come from step 4. We assume sy = 0.17 [56] independent of 
X. Parameters sz and pyz are the unknowns to be determined. With the first guess 
of Sz from step 4 and an initial guess of pyz (e.g. 0.4) we generate our ffist mock set 
including the relevant correlation. This mock set reproduces the SMF because we are using 
the bias corrected abundance matching X-Y relation. However, the set does not reproduce 
the observed VDF and the observed stellar mass-velocity dispersion relation because of the 
inaccurate sz and the unknown pyz- 

6. We use the mock set itself to estimate the bias in sz- This is done in the following 
way. First, calculate the standard deviation of Z at fixed Y, sz{Y). Next, take the lower 
dimensional set {X, Y} from the mock set and then assign Z toY using sz{Y) just calculated, 
pretending the correlation in the mock set were unknown. Thus, we get a biased set {X, Y, Z}. 
Finally, we calculate sz at fixed X from this biased set and then compare it with that of the 
mock set to estimate the bias. 

7. Correct sz for the bias estimated in step 6. Using this corrected sz generate a revised 
mock set {X, Y, Z} using the bivariate normal probability distribution of Y and Z at fixed 
X. We can check that this revised mock set results in an improved match to the observed 
VDF. We note that the resulting VDF has much to do with sz but little to do with pyz- The 
correlation pyz has much to do with the resulting Y-Z relation, its mean relation and the 
dispersion of Z at fixed Y. Hence, for the given sz{X) we adjust pyz so that the resulting 
Y-Z relation matches the observed stellar mass- velocity dispersion as closely as possible. For 
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Figure 6. Parameters of the bivariate normal distribution of Y [= logiQ{Mi, / Mq)] and Z 
logig(CT/km s~^)] as functions of X [= \ogiQ{M^i,- /Mq)]. 



the sake of simplicity we allow pyz to be only a linear function of X with broken slopes. 

8. If the results of step 7 are not satisfactory go back to step 6 and iterate. Initially 
we focus on reproducing the VDF by correcting sz{X) with a constant pyz- Eventually the 
initially adopted pyz will be modified as the observed median relation between stellar mass 
and velocity dispersion cannot be well reproduced with a constant pyz- We find that with a 
constant pyz = 0.35 the resulting VDF matches well the observed VDF for a > 100 km s"^. 
This means that through the iteration we have found a proper functional behavior of sz{X) 
(see Fig. 6). On the other hand, the resulting mean value of Z at Y deviates downward 
from the observed relation towards the high mass end. This requires us to correct the 
constant pyz- After trial and error we find that a varying pyz{X) with enhanced values 
for relatively more massive systems, as shown in Fig. 6, can reproduce the observed median 
relation between stellar mass and velocity dispersion (see Fig. 5). However, for this model 
of pyz{X) the intrinsic scatter of Z at fixed Y is slightly underestimated as can be seen 
in Fig. 5. This lowered dispersion is a consequence of the strengthened correlation. This 
manifests the difficulty of matching the VDF and the stellar mass-velocity dispersion relation 
simultaneously to a high precision. Some compromise is necessary and we have chosen to give 
more weight to the VDF for a > 100 km s~^ and the median stellar mass-velocity relation 
but less weight to the the intrinsic scatter of velocity dispersion at fixed stellar mass as the 
observational error of the latter has not been quantified. 
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3.2 The result 

We obtain a final set {Mvir, M*, a} of ~ 6.9 x 10® mock early-type systems with Myir > 
IO^^'^Mq for a training comoving volume of 4 x 10^ Mpc^. A subset of 50,000 systems 
is displayed in Fig. 7. The projected Mvir-M^ relation is compared with recent available 
observational results based on other methods including galaxy-galaxy weak lensing [59] , halo 
occupation statistics [42], and kinematics of satellite galaxies [56]. Our mean Myir-M^ relation 
lies near the median of other results for Mvir < W^^Mq but lies lower down to —0.2 dex for 
Mvir > 10^^ Mq. We note that this high mass part behavior of the abundance matching Mvir- 
M* relation in comparison to other methods is also found for all-type galaxies [49] . However, 
it turns out that this difference has relatively minor effect on our results on dark matter 
density profiles. The Mvir-o" relation for early-type galaxies is derived here for the first time. 
Abundance matching Mvir-o" relations for all- type galaxies can be found in [50, 60]. 



4 Jeans dynamical modeling of early-type galaxy-halo systems 

Through the bivariate normal distribution of log]^Q(Mv,) and log]^o(^) a function of Mvir 
for early-type spheroidal galaxy-halo systems, stellar mass M^, and stellar velocity dispersion 
a are assigned to the galaxy embedded at the the center of each halo with virial mass Mvir. 
The galaxy with M^ and a assigned can be further assigned its stellar mass density profile 
through the observed correlations with the effective radius Re and the Sersic index n. These 
correlations are described in appendix A. 

4.1 The spherical Jeans equation 

Once the velocity dispersion and the parameters of the stellar mass distribution are specified, 
we can use the Jeans equation to constrain the dark matter density profile of the embedding 
halo because the velocity dispersion depends not only on the stellar mass distribution but 
also on the dark matter distribution. We focus on the spherical Jeans equation [61] given by 

«^ + 2^[Mr)aHr)] = -G^^^^^f^^^ (4.1) 
dr r 

where C7r(r) is the radial stellar velocity dispersion at radius r, /9*(r) is the three-dimensional 
stellar mass density at r, M(r) is the total (i.e. stellar plus dark) mass enclosed within r, i.e. 
M(r) = M*(r) -|- Mdm(?'), and /3(r) is the velocity dispersion anisotropy at r given by 

where (7g{r) and cr^{r) are the tangential velocity dispersions in spherical coordinates. An 
integral solution of the Jeans equation for crrif) is given in appendix B. 

The LOSVD of stars at projected radius R on the sky (Tios(i?) is given by 



oo 



R2 



2 



1 - 



Vr2 - i?2 ' 



(4.3) 



where S^(i?) is the two-dimensional stellar mass density projected on the sky. The stellar 
mass weighted j-th power of LOSVD within an aperture of radius R is then given by 
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Figure 7. (Top) The bivariate normal distribution of log]^o(Ai^/M0) and \o%^^^{a jVva. s~^) as a 
function of log]^Q(Mvir/-/W0) for early-type galaxy-halo systems where Mvir is the virial mass of the halo 
embedding an early-type galaxy at the center and M-i, and a are the stellar mass and central stellar 
velocity dispersion of the galaxy respectively. (Bottom) Two projections of the three-dimensional 
distribution. Black curves are the mean relations. The mean Mvir-M^ relation is compared with 
recent observational results on early-type or red galaxies based on galaxy-galaxy weak lensing (blue 
curve) [59], halo occupation statistics (green) [42] and satellite kinematics (cyan) [56]. 



Finally, assuming that stellar mass follows luminosity the quantity to match the SDSS velocity 
dispersion a is given by 



Ks)(ii = iic/8). 



(4.5) 
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4.2 Models for the dark matter distribution 



Initially all mass is supposed to follow a density profile similar to the NFW profile [62-64], 
as predicted by N-body simulations. As the galaxy is formed and settled at the center of 
the halo, the separate stellar mass distribution is embedded in the remaining dark matter 
distribution. While the stellar mass distribution can be represented by the Sersic profile as 
indicated by observations, the remaining dark matter distribution is completely unknown. 

With the assumption that the dark matter distribution has been smoothly readjusted 
from the initial NFW(-like) profile, we model the unknown dark matter distribution using a 
generalized NFW density profile, referred to as aNFW throughout, given by 



Ps 

{r/nni + r/nY 



where a is the inner density power-law slope (a = 1 being the NFW value), rg is the scale 
radius related to the concentration parameter Cvir via 

Cir = ^, (4.7) 

and the parameter is related to the dark matter mass within the virial radius -/Vfdm(= 
Mvir — M^) via 

^^-4^r|/„(c„,)' ^^-^^ 

where the function fa{x) is given by 

We also consider the Einasto profile [65], that is suggested by relatively more recent dissipa- 
tionless N-body simulations [63, 64], given by 

PEinir) = p^2 exp {-(2/a) [(r/r„2)° - l] } , (4.10) 

where r_2 is the radius at which the logarithmic slope of the density is —2 and from which 
we define another concentration parameter given by 



(4.11) 



and the parameter /9_2 is related to the dark matter mass via 

4vrr32/a(cvir) 



P-'i = Z ? ,~ X ' i^-^^) 



where the function fa{x) is given by 



h{x) = 4 ( lY'\^V 1,4^" ) (4.13) 



a \ a I \ <y J \ a a 



where 7(0, b) is the incomplete gamma function. 
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Figure 8. (Left) Constraints on the aNFW halo dark matter density profile parameters a and 
Cvir with either of them fixed at the Bolshoi N-body simulation predicted NFW value based on the 
spherical Jeans equations with constant anisotropy of velocity dispersions. Solid curves are the median 
while dashed curves are the 68% probability limits. The green curves are based on increased stellar 
masses by a factor of 2 for Mvir > IO^'^A/q and a factor from 1 to 2 for lO^^Afg > Mvir > IO^Mq. 
(Right) Upper panel shows the simultaneously constrained values of constant anisotropy /3 with the 
prior [—0.5,0.5]. Lower panel shows the minus logarithmic slope of the predicted total (i.e. stellar 
plus dark matter) mass density profile at r = i?e/2. 



4.3 Results on dark matter density profiles 

4.3.1 Restricted results based only on the SDSS velocity dispersion 

The dark matter distribution given by the aNFW (or Einasto) profile is completely specified 
by two parameters a and Cvir (or, a and Cvir)- With the stellar velocity dispersion and the 
stellar mass distribution assigned, the Jeans equation constrains the plane spanned by a and 
Cvir provided the velocity dispersion anisotropy /3(r) is specified in advance (or constrained 
simultaneously in conjunction with other constraints). This is done by equating the theo- 
retical velocity dispersion given by equation (4.5) to the SDSS velocity dispersion. Let us 
consider the case of constant anisotropy and assign its value using uniform random variable 
from the range [—0.5,0.5] [6, 7, 10, 66]. Let us further consider the case that either a or Cvir 
is fixed so that the other parameter is determined by the Jeans equation for a given value of 

Fig. 8 shows the results for the aNFW profile as a function of the halo mass M^-„. For 
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the first case Cvir is fixed at the value predicted by the Bolshoi N-body simulation. In this case 
the constrained values of a lie above the NFW value a = 1 with increasingly larger values as 
Mvir decreases. For the second case, a is fixed at the NFW value and Cvir is constrained by the 
Jeans equation. Similarly to the first case the values of Cvir lie above the N-body simulation 
prediction. These results indicate that dark matter densities are boosted in the inner halos of 
early-type galaxies. The density boost is greater for a lower mass halo. Qualitatively similar 
results are obtained for the Einasto profile as well. If a is fixed at 0.17 (the N-body predicted 
value) [63], the constrained values of Cvir tend to lie above the Bolshoi prediction. If Cvir is 
fixed at the value predicted by the Bolshoi simulation, the constrained values of a tend to 
lie below a = 0.17 implying steepened inner density slopes. 

These results are for specific stellar masses of the embedded galaxies that are based on 
the assumption of the Chabrier initial mass function (IMF) [41]. Stellar masses would be 
larger by a factor of two (a dex of 0.3) if the Salpeter IMF were assumed instead. Furthermore, 
as shown in Fig. 7 our statistical match between galaxies and halos gives a lower by up 
to 0.2 dex at fixed Myir for Mvir > IO^^Mq compared with other independent results that 
are based on IMFs similar to the Chabrier. To gauge the possible systematic effect due to 
stellar mass uncertainties we consider adjusting stellar masses by 0.3 dex for Myir > 10^^ Mq 
and O.3[logio(Mvir/M0) - 13] for IO^^Mq > M^h- > IO^^Mq. For the case of using the fixed 
NFW inner slope a = 1 Cvir is lowered by up to ~ 30 percent. 

The above two cases (fixing a/a or Cvir/cvir) imply different inner density profiles that 
cannot be distinguished by the Jeans equation alone. It shows the degeneracy of the allowed 
density profile at a given Mvir. It can be seen by the bottom right panel of Fig. 8 which 
shows the constrained values of the negative logarithmic slope 7 of the three-dimensional total 
radial density profile at Rc/2, i.e. 7 = — (iln[p^(r) + P(i^{r)]/dlnr. In fact, the degeneracy 
is broader than the above two cases. In other words, other combinations of a and Cvir (or, a 
and Cvir) can satisfy the Jeans equation as well. To break the degeneracy we need additional 
constraints. Those constraints may also help to constrain the velocity dispersion anisotropy. 

4.3.2 Dark matter density profiles constrained by velocity dispersion profiles 

The SDSS velocity dispersion corresponds to the luminosity- weighted LOSVD within Rc/8. 
Clearly, this alone cannot constrain uniquely the dark matter distribution in the halo. Mea- 
sured velocity dispersions at multiple radii for each SDSS galaxy would be most useful in 
breaking the degeneracy. Without them we consider independent observational constraints 
on the velocity dispersion profile (VP) of early-type galaxies [8, 67]. The available obser- 
vational constraints show that the logarithmic slope r] of the luminosity-weighted LOSVD 
profile at Rc/2, i.e. ?? (iln(crios)(^)/c?ln-R, ranges from -0.066 ± 0.035 [8] to -0.04 [67]. 
We assume that the mean value of r] is given by the mean of these two independent mea- 
surements, i.e., {rj) = —0.053. For the dispersion of r] we assume the root-mean-square of the 
measured dispersion 0.035 for one sample and the difference between the two independent 
means 0.026, i.e., ar, = 0.044. 

We then proceed as follows. For a halo of Mvir with a and stellar mass distribution 
assigned, we first obtain a constrained set of {a,Cvir,/3} using equation (4.5) only. We then 
use the above VP constraints to select the most likely model out of the degenerate set. We 
select our best-fit model by minimizing the figure-of-merit function defined by 




(4.14) 
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Figure 9. Constraints on the same parameters of Fig. 8 with a and Cvir simuUaneously aUowed to 
vary based on velocity dispersion profile constraints as well as the spherical Jeans equation. The blue 
solid curves show the expectation values. Orange, yellow and gray regions contain the 68%, 95% and 
99.7% occurrences respectively. The green solid and dashed curves show the expectation value and 
the 95% limits based on increased stellar masses as in Fig. 8. The black dashed line in the lower left 
panel shows the mean NFW concentration predicted by the Bolshoi N-body simulation. The black 
dashed line in the lower right panel corresponds to the isothermal profile. 



where ■rf^°^'^^ is the VP slope of each model system at Re/2. If the value of Q^yp) ^® greater 
than 9 for the best-fit model, then the whole degenerate set is rejected. But, this occurs 
rarely. 

Fig. 9 shows the results on the distribution of the dark matter density profiles based on 
the aNFW model with constant velocity dispersion anisotropy in the range [—0.5,0.5]. The 
constrained ranges of a and Cvir show that their expectation values are decreasing functions 
of Mvir but there are significant halo-to-halo variations. For the majority of occurrences 
with Mvir < 10^^ Mq the values of a lie above the NFW value. For 68% of galactic halos 
with lO-^^M© < Mvir < lO^^-^M© we have 1 < a < 1.7. The values of Cvir scatter around 
predictions of cosmological dark matter simulations [44, 53]. 

Fig. 9 also shows the distribution of the negative logarithmic slope 7 at Rc/2. For 
halos with 10^^ Mq < Mvir < 10^^ Mq the values of 7 scatter around the isothermal value 
7 = 2 with the expectation value remarkably close to the isothermal. This result for a 
general sample of early-type galaxies is in good agreement with the results for strong lensing 
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Figure 10. Same as Fig. 9 except that the Einasto model parameters dt and Cvir are constrained. 
Parameter 7dm is the minus logarithmic slope for the Einasto model dark matter distribution, i.e. 
7dm ^ —dh\ p^ii^{r) / dlnr at Rc/2. The constrained values of 7dm match well with those of a of the 
aNFW model shown in Fig. 9. 



observations of early-type galaxies [9-11]. For massive cluster-sized halos the central total 
density profile within the optical region tends to be shallower than the isothermal. The 
distribution of the velocity dispersion anisotropy f3 is tilted toward a mild radial anisotropy 
with the expectation value close to isotropy. This distribution is consistent with independent 
results for individual systems [6, 10, 66]. 

Fig. 10 shows the corresponding results for the Einasto model. The constraints on f3 and 
7 are strikingly similar to those for the aNFW model. The constraints on a and Cvir cannot be 
directly compared with those on a and Cvir- For the Einasto model the logarithmic slope of the 
density varies continuously toward the origin, i.e. 7dm = 2(cvirr/rvir)" for Pdm('") oc r~'^'^'^, 
so that unlike the aNFW model an inner density slope is not characterized by a single 
parameter. At a fixed inner radius the slope 7dm can be boosted either by increasing Cvir 
or decreasing d. Fig. 10 indicates that the former is likely the case for galactic halos. The 
constrained values of 7dm at -Re/2 are strikingly similar to those of a of the qNFW model. 
A remarkable feature of Fig. 10 is that the constrained values of d center around the N-body 
predicted value of dNbody = 0.17 insensitive to halo mass, despite the result that the dark 
matter density slope is boosted in the inner optical region. 

We have tested our results against a couple of dynamical mass scaling relations presented 
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recently in the literature. One scaling involves the mass within r ~ i?e/2 and the LOSVD 
at that radius [68] while the other scaling involves the mass within r ~ 4i2e/3 and the 
luminosity- weighted LOSVD to infinity [69]. Prom our results the average masses predicted 
by these scalings are accurate within 10% for systems with lO""^^ Mq < Mvir < 10"^^ Mq and 
there are system-to-system variations with a typical dispersion of 5% to 10% at fixed Mvir 
(appendix C). 

5 Implication for bar yon- induced halo contraction and characterization of 
the contracted profiles 

Our constrained dark matter density profiles of the halos support the hypothesis that halos 
are modified in response to dissipational gas cooling and galaxy formation [39, 40]. Our results 
indicate that most halos have contracted and dark matter density profiles of the contracted 
halos are diverse. Fig. 11 shows the ratios of the dark matter masses of our halos to those of 
NFW halos with the average concentrations predicted by the Bolshoi dissipationless N-body 
simulation within three radii, i.e. the effective radius Re, the NFW scale radius rs,NFW and one 
fifth of the halo virial radius (rvir/5). If the NFW profile were assumed for the initial N-body 
predicted halo and there were no modification to the halo due to baryonic effects, then these 
ratios should center around unity with scatters arising solely from the scatters of N-body 
predicted NFW concentrations. Fig. 11 shows evidently dark matter density enhancement 
in the optical region of the inner halo. On average, the density enhancement factor within 
r = i?e increases steeply from ~ 1 (no enhancement) for Mvir > 10^^ Mq to 4 - 5 for 
Mvir < 10^^ Mq. For lower mass halos, however, the scatter also increases rapidly indicating 
greater diversity. At the scale of Ts^nfw (i-e- the N-body predicted scale radius assuming the 
NFW model) the halo density enhancement is weakly present only for Mvir < lO^'^'^ Mq. At 
the scale of rvir/5 there is no longer statistically significant density enhancement. A statistical 
analysis of weak lensing effects for r > rvir/5 based on statistically representative samples of 
galaxies, groups and clusters from SDSS yields concentrations similar to N-body predictions 
assuming the NFW profile for lO^^j^^ < j^^.^ < lO^^M© [70]. This result is consistent with 
our finding of no or little density enhancement within r = rvir/5. The weak lensing measured 
scatter of concentration gives rise to Icr scatter of 10% - 20% in mass within r = rvir/5 
consistent with our result shown in Fig. 11. 

Fig. 11 also shows the fraction of dark matter mass within a sphere of radius r = for 
our halos embedding early-type galaxies. It is weakly dependent on Mvir. For the majority 
of our constrained systems dark matter fraction within R^ lies between 20% to 70%. Dark 
matter fraction at fixed Mvir covers a broad range showing a great diversity. Dynamical 
analyses of nearby elliptical galaxies indicate ~ 30% of dark matter contribution within Rg 
[66] and a recent analysis of strong lensing galaxies indicates that fdm{'>' < Rc) can be ~ 0.2 
- 0.6 [11]. These results agree well with our results. 

Traditionally, adiabatic contraction has been invoked for the response of halos to dis- 
sipational baryonic processes [24, 25, 39, 40]. Hydrodynamic simulations have been used to 
study adiabatic contraction in detail recently [24, 25]. Independent of hydrodynamic sim- 
ulations and relying only on SDSS observed galaxy data and robust classical results from 
N-body dissipationless simulations, our results demonstrate evidently that halo contraction 
takes place in most occurrences of spheroidal galaxy formation. 

An observational indication of adiabatic contraction in a statistically representative 
sample of SDSS elliptical galaxies is reported recently [20] based on a combination of weak 
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Figure 11. (Left panels & Top Right) Ratio of the dark matter masses within a certain three 
dimensional radius for the initial NFW halo without galaxies (Mnfw) and the constrained realistic 
halo embedding an early-type galaxy at the center (Mdm) based on the aNFW model. Here Ts^nfw 
is the scale radius of the NFW profile while is the effective radius of the embedded central galaxy. 
These results show that most halos have contracted in the optical region of central galaxies with the 
varying degree of mass enhancement depending on the the halo virial mass. However, the contraction 
has modified the halo up to the scale radius only for relatively less massive halos with Mvir < lO^'^M© 
and starts to be unimportant beyond the scale of 20% of the virial radius. (Bottom Right) Fraction of 
dark matter mass within i?c of the early-type galaxy embedded at the center of a halo. Green curves 
show the results based on the increased stellar masses as in Fig. 9. The corresponding results for the 
Einasto model are similar to those shown here. 



lensing measurements of halos and stellar velocity dispersions. These authors assume the 
NFW profile and notice that their analysis requires higher concentrations of halos than pure 
dark matter simulation halos. A combined analysis of strong and weak lensing effects for 28 
clusters from the Sloan Giant Arcs Survey finds increasingly over-concentration as halo mass 
gets lower assuming the NFW profile [23]. 

Based on generalized dark matter density profiles that are allowed to vary from the NFW 
profile, our results provide not only evidence for halo contraction but also a characterization 
of contraction. Our results indicate that contracted halo profiles generally deviate from the 
NFW and display a range of diversity and halo contraction is more frequently realized by 
steepening of the radial density slope rather than rescaling of the NFW profile. If the NFW 
profile were forcibly used to describe the contracted halo, then the concentration would 
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be significantly higher than that of the pure dark matter simulation halo which explain 
high concentrations of halos obtained by recent observational studies using the NFW profile 
[20, 23]. 



6 Characterizing dark matter annihilation strength in the halos of early- 
type galaxies 

One proposed way of identifying dark matter particles has been the detection of 7-ray emis- 
sion from dark matter annihilation or decay in the halos [19, 71]. In the case of dark matter 
annihilation the predicted 7-ray flux is known to be very sensitive to dark matter density. 
Specifically, the predicted 7-ray flux from dark matter annihilation is proportional to the dark 
matter density squared integrated along the line-of-sight. When averaged over the angular 
scale of the halo it can be approximated by [72] 



where D is the distance to the halo. 

Clearly, halo contraction is expected to boost the predicted 7-ray flux from dark matter 
annihilation. Fig. 12 shows the predicted values of Jdm x for our constrained halo dark 
matter distributions based on the Einasto model. The aNFW model cannot be used for 
calculating Jdm because Jdm diverges for a > 1.5 as r — )• 0. The specific values of Jdm 
at a fiducial distance oi D = d^A Mpc are also shown and compared against several local 
dwarf spheroidal galaxies [73] and low redshift clusters [72] that have been proposed as 
promising targets for 7-ray flux. Fig. 12 suggests that typical nearby early-type galaxies with 

can also be promising targets. This rather surprising prediction 
is the consequence of the mass dependent enhancement of Jdm ^is shown in the middle panel 
of Fig. 12, resulting from the varying degree of halo contraction. We note that there are two 
well-known elliptical or lenticular galaxies at a distance of ~ 4 Mpc, Centaurus A (NGC 5128) 
and Maffei 1. Although our results are only for the halos embedding early- type galaxies at 
their centers, they highlight the importance of halo contraction induced by galaxy formation 
for dark matter search. 

7 Conclusions 

Based on a combination of the observed statistical properties of SDSS galaxies and the halo 
mass function from the Bolshoi N-body dissipationless simulation in conjunction with the 
observed properties of velocity dispersion profiles of early-type galaxies, we find for the halos 
embedding early-type galaxies: 

• The inferred dark matter density profiles of the halos imply significantly higher dark 
matter densities in the inner regions compared with those from the dissipationless 
simulation. The mean density boost factor ranges from ^ 1 for M^ir > lO^^M© to 
~ 4 — 5 for Mvir < IO^^Mq. This provides an independent evidence for halo contraction 
for galactic halos. 

• Using a couple of general-class models for the dark matter distribution we obtain sta- 
tistical characterizations of the contracted profiles. The three-dimensional dark mat- 
ter density slope 7dm at r = Rc/2 for /5dm(^) oc r~'^'^'^ is found to be in the range 
1 < 7dm < 1-7 for 68% of galactic halos with lO^^jvi^ < j^^^.^ < IO^^-^Mq. 




(6.1) 
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• Halo contraction boosts significantly dark matter annihilation strength in the halos 
embedding early-type galaxies so that nearby early-type galaxies may be promising 
targets for indirect dark matter search. 
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A SDSS early-type galaxy parameter correlations 

An early-type galaxy with and a specified may be given its stellar mass density profile 
using observed correlations with R^. (effective radius) and n (Sersic index) assuming the Sersic 
stellar mass density profile [54] given by 



where X = R/Re {R being the two dimensional radius), 5^ = 2n — 1/3 + 0.009876/n (0.5 < 
n < 10), and An = b'^ / [2TTnT {2n)] [74, 75]. This is the projected two dimensional density 
normalized such that its integrated total mass is unity. The deprojected three dimensional 
density at x = r/R^ (the three dimensional radius in units of Rq) is then given by 



with Q„ = 1 - 0.6097/n + 0.05563/n2 and />„ = (6„)"(3-«n)/{47rnr[n(3 - a„)]} [74, 75]. 



To a galaxy with and a specified we give a specific value of Rc using the observed 
distribution of M^y^/M^ as a function of M^, where -/Vf^yn = bR^a'^/G {G being Newton's 
gravitational constant) is a 'dynamical' mass. Fig. 13 shows the observed distribution of 
M(jyn/M^. Notice that the distribution of Mjyn/M^ essentially shows the scatter of the 
fundamental mass plane relation of early- type galaxies for three parameters R^, cr and M^. 

Finally, we give a specific value of n to a galaxy using the observed distribution of n as a 
function of (see Fig. 14). This result comes also from SDSS data [42]. Sersic index n may 
also be correlated with a and Rq. However, the current unavailability of such correlations 
forces us to use only the correlation shown in Fig. 14. 

B An integral solution of the Spherical Jeans equation 

For a specific form of the velocity dispersion anisotropy /3(r) the radial velocity dispersion 
crr(r) can be expressed as an integral form that can be easily evaluated numerically. 




(A.l) 




(A.2) 
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Figure 12. (Top) The dark matter density squared integrated along the line-of-sight and averaged on 
the halo angular scale, Jdm times the halo distance squared D^, approximated by JdmD^ — J p'^{r)dr 
[72]. The results are shown only for the Einasto model because Jdm diverges for certain central 
cusps of the aNFW model. The predicted 7-ray flux from dark matter annihilation is proportional 
to Jdm- (Middle) Ratio of Jdm for the constrained halo embedding a central early-type galaxy to 
that for the dark matter only N-body simulation halo with a = 0.17 and Cvir matched to the Bolshoi 
concentrations. This shows that the predicted 7-ray flux is enhanced in most galactic halos embedding 
spheroidal galaxies. (Bottom) The predicted range of Jdm for nearby normal spheroidals (i.e. ellipticals 
and lenticulars) at a fiducial distance oi D = d^A Mpc is shown and compared with the estimates for 
local dwarf spheroidal galaxies (blue arrows) [73] and low-redshift (2 < 0.025) clusters (red arrows) 
[72]. 



We recast the spherical Jeans equation [eq. (4.1)] as 

dy{r 



dr 



+ P{r)y{r) = q{r), (B.l) 
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Figure 13. (Top) The distribution of logiQ(Mcjyn/A/*) for SDSS early-type galaxies with Afdyn = 
5Re(T^/G. We obtain this distribution from the data analyzed in [41]. Data points are the median 
values and 68% scatters. Red solid curve and dashed curves are the polynomial fits assuming the 
Gaussian distribution. (Bottom) The distribution of the mean value of log]^g(i?e/kpc) based on the 
Gaussian mean of logjo(-^^dyn/Af*) from the top panel and the mean of log]^Q(cr/km s"-'^) from Fig. 5. 



where y(r) = pi,{r)a^{r), p{r) = 2/5(r)/r, and q{r) = —Gpi,{r)M{r)/r'^. Then, with the 
definition of Lo{r) = exp ^J^ p{t)dt^ we get 



y(r) 



uj{r) 



oj{t)q{t)dt. 



(B.2) 



Let us now consider a general functional form for /3(r) given by [69] 



y.2 _|_ j,2 



(B.3) 
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Figure 14. The distribution of Sersic index (n) as a function of stellar mass (M-^) for SDSS early- 
type galaxies (red points) with the curves showing median values (full) and 68% scatters (dashed). 
These data are taken from [42]. 



where /3o and /3oo correspond respectively to the anisotropies at the origin and infinity and 
is the transition radius between the two. For this model uj{r) becomes 

uj{r) = r2*(r2 + r^Y^o^-M^ (B.4) 

Substituting this into equation (B.2) we get 

\ J^2^/3o— /3oo />oo 

a^r) = t/,, / t'^'^-'^i^' + r}r--'°^Mt)M{t)dt. (B.5) 
For the case of constant anisotropy /3 = /3o = (3oo it takes the following simple form 

-2/3 foo 

a^{r) = G—- / t^^^~^^p^{t)M{t)dt. (B.6) 

In the above the total mass within r is given by M(r) = M^(r) + M^-^{r) with the 
stellar mass 

^ ,, j[«(3- . ),i,.(r//4)-/"] 

r[n(3 - a„)J 
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and the dark matter mass 



Md^(r) = (M.i, - M^-f^ (B.8) 

for the aNFW model and 

Md^(r) = (M^ir - M,) ^fl^~^^ , (B.9) 

for the Einasto model. In the above T{x) and 7(x, y) are the gamma function and the incom- 
plete gamma function respectively. The functions fa{x) and fa{x) are given by equation (4.9) 
and equation (4.13) respectively. 



C Dynamical mass scaling relations and alternative results for halo mass 
profile 

Our results on the halo profiles presented in the main text are only based on the observa- 
tional constraints on the velocity dispersion profile (VP). Dynamical mass scaling relations 
have been presented recently in the literature [68, 69]. These relations are derived from the 
spherical Jeans equation in conjunction with physically well-motivated dynamical assump- 
tions that are consistent with current observations. These relations allow us to estimate 
dynamical masses within certain radii from velocity dispersions. A crucial common feature 
of these relations is that they are insensitive to the value of the velocity dispersion anisotropy 
p. 

The first mass estimate [68] is expressed as 

M(-*)(ropt) = ^roptaL(ropt), (C.l) 

where ropt is an optimal radius defined in the reference and ropt ~ 0.5i?c typically. The 
numerical factor /„ is a velocity ratio defined in the reference and typically ~ 0.6. This 
mass estimate is based on the assumption that the total gravitational potential is that of 
the isothermal potential. This is not strictly satisfied by observed galaxies but a reason- 
able approximation for the potential within stellar extents of galaxies supported by various 
astrophysical studies including strong lensing [9-11]. 
The second mass estimate [69] is expressed as 

M(°^*Hri/2) = |ri/2Ks)(00), (C.2) 

where ri/2 is the three-dimensional radius at which the enclosed stellar mass is and 
~ 4i2e/3 for most stellar mass distributions of spheroids. This mass estimate is based on 
the assumption that (j\os{R) varies sufficiently slowly with R near i?e as supported by current 
observations [66, 67]. 

Because of the approximate assumptions made we expect that our galaxy-halo systems 
will satisfy the above mass scaling relations only approximately. Fig. 15 shows the ratios of the 
estimated masses based on the scaling relations to the masses enclosed within the relevant 
radii for our systems. For lO^^ < m^-,^ < 10^^ M© the mean estimated masses agree 
with the mean enclosed masses within 10%. However, the mean ratio (M^^*^*^ (ropt )/-/W (ropt)) 



- 26 - 




logio(MviyM^) 

Figure 15. Ratios of the estimated masses based on two scaling relations from the literature to the 
masses within two radii rj/2 (~ 4_Ro/3, left panel) and ropt -Rc/2, right panel) for our galaxy-halo 
systems constrained by the observational VP constraints. The orange, yellow and gray regions include 
the 68%, 95% and 99.7% of our systems respectively. The blue curves are the average values. 



(blue curve in the left panel of Fig. 15) tends to lie above unity whereas the mean ratio 
(M'-'^^*)(ri/2)/-^('"i/2)) (blue curve in the right panel of Fig. 15) systematically varies as a 
function of M^i^-. We do not speculate on possible sources of these behaviors. 

In the above we have tested the dynamical mass scaling relations with our systems 
constrained by the VP constraints only. Alternatively, we could use the dynamical mass 
scaling relations as additional constraints on the halo dark matter density profiles. We 
proceed by introducing two figure-of-merit functions defined by 

and 




(WoU) - j2 I M{r,/2) 



where we assume 6 = 0.1. We then combine these functions to the function Qjyp^ based 
on the VP constraints given by equation (4.14). Fig. 16 and Fig. 17 show the results on 
the aNFW model based on combined figure-of-merit functions given by Q^yp) + Q(chur) 
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Figure 16. Same as Fig. 9 except that this result is based on the combined figure-of-mcrit function 
of Q(vp) +'3(Chur) [equation (4.14) & equation (C.3)]. 



and Q(vp) ~^ Qfwoii) • These results are not particularly tighter or statistically significantly 
different from the result based on the VP constraints only shown in Fig. 9. 
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